Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Fine-mapping genomic loci refines bipolar disorder risk genes.

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Clinical, genomic, and neurophysiological correlates of lifetime suicide attempts among individuals with alcohol dependence.

Research has identified clinical, genomic, and neurophysiological markers associated with suicide attempts (SA) among individuals with psychiatric illness. However, there is limited research among those with an alcohol use disorder (AUD), despite their disproportionately higher rates of SA. We examined lifetime SA in 4,068 individuals with DSM-IV alcohol dependence from the Collaborative Study on the Genetics of Alcoholism (23% lifetime suicide attempt; 53% female; 17% Admixed African American ancestries; mean age: 38). We 1) conducted a genome-wide association study (GWAS) of SA and performed downstream analyses to determine whether we could identify specific biological pathways of risk, and 2) explored risk in aggregate across other clinical conditions, polygenic scores (PGS) for comorbid psychiatric problems, and neurocognitive functioning between those with AD who have and have not reported a lifetime suicide attempt. The GWAS and downstream analyses did not produce any significant associations. Participants with an AUD who had attempted suicide had greater rates of trauma exposure, major depressive disorder, post-traumatic stress disorder, and other substance use disorders compared to those who had not attempted suicide. Polygenic scores for suicide attempt, depression, and PTSD were associated with reporting a suicide attempt (ORs = 1.22-1.44). Participants who reported a SA also had decreased right hemispheric frontal-parietal theta and decreased interhemispheric temporal-parietal alpha electroencephalogram resting-state coherences relative to those who did not, but differences were small. Overall, individuals with alcohol dependence who report SA appear to experience a variety of severe comorbidities and elevated polygenic risk for SA. Our results demonstrate the need to further investigate suicide attempts in the presence of substance use disorders.

Genetic architecture distinguishes tinnitus from hearing loss.

Tinnitus is a heritable, highly prevalent auditory disorder treated by multiple medical specialties. Previous GWAS indicated high genetic correlations between tinnitus and hearing loss, with little indication of differentiating signals. We present a GWAS meta-analysis, triple previous sample sizes, and expand to non-European ancestries. GWAS in 596,905 Million Veteran Program subjects identified 39 tinnitus loci, and identified genes related to neuronal synapses and cochlear structural support. Applying state-of-the-art analytic tools, we confirm a large number of shared variants, but also a distinct genetic architecture of tinnitus, with higher polygenicity and large proportion of variants not shared with hearing difficulty. Tissue-expression analysis for tinnitus infers broad enrichment across most brain tissues, in contrast to hearing difficulty. Finally, tinnitus is not only correlated with hearing loss, but also with a spectrum of psychiatric disorders, providing potential new avenues for treatment. This study establishes tinnitus as a distinct disorder separate from hearing difficulties.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders.

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Single-Nucleus Transcriptome Profiling of Dorsolateral Prefrontal Cortex: Mechanistic Roles for Neuronal Gene Expression, Including the 17q21.31 Locus, in PTSD Stress Response.

OBJECTIVE: Multidisciplinary studies of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) implicate the dorsolateral prefrontal cortex (DLPFC) in disease risk and pathophysiology. Postmortem brain studies have relied on bulk-tissue RNA sequencing (RNA-seq), but single-cell RNA-seq is needed to dissect cell-type-specific mechanisms. The authors conducted the first single-nucleus RNA-seq postmortem brain study in PTSD to elucidate disease transcriptomic pathology with cell-type-specific resolution. METHOD: Profiling of 32 DLPFC samples from 11 individuals with PTSD, 10 with MDD, and 11 control subjects was conducted (∼415K nuclei; >13K cells per sample). A replication sample included 15 DLPFC samples (∼160K nuclei; >11K cells per sample). RESULTS: Differential gene expression analyses identified significant single-nucleus RNA-seq differentially expressed genes (snDEGs) in excitatory (EX) and inhibitory (IN) neurons and astrocytes, but not in other cell types or bulk tissue. MDD samples had more false discovery rate-corrected significant snDEGs, and PTSD samples had a greater replication rate. In EX and IN neurons, biological pathways that were differentially enriched in PTSD compared with MDD included glucocorticoid signaling. Furthermore, glucocorticoid signaling in induced pluripotent stem cell (iPSC)-derived cortical neurons demonstrated greater relevance in PTSD and opposite direction of regulation compared with MDD, especially in EX neurons. Many snDEGs were from the 17q21.31 locus and are particularly interesting given causal roles in disease pathogenesis and DLPFC-based neuroimaging (PTSD: ARL17B, LINC02210-CRHR1, and LRRC37A2; MDD: LRRC37A and LRP4), while others were regulated by glucocorticoids in iPSC-derived neurons (PTSD: SLC16A6, TAF1C; MDD: CDH3). CONCLUSIONS: The study findings point to cell-type-specific mechanisms of brain stress response in PTSD and MDD, highlighting the importance of examining cell-type-specific gene expression and indicating promising novel biomarkers and therapeutic targets.

Differences in genetic correlations between posttraumatic stress disorder and alcohol-related problems phenotypes compared to alcohol consumption-related phenotypes.

BACKGROUND: Posttraumatic Stress Disorder (PTSD) tends to co-occur with greater alcohol consumption as well as alcohol use disorder (AUD). However, it is unknown whether the same etiologic factors that underlie PTSD-alcohol-related problems comorbidity also contribute to PTSD- alcohol consumption. METHODS: We used summary statistics from large-scale genome-wide association studies (GWAS) of European-ancestry (EA) and African-ancestry (AA) participants to estimate genetic correlations between PTSD and a range of alcohol consumption-related and alcohol-related problems phenotypes. RESULTS: In EAs, there were positive genetic correlations between PTSD phenotypes and alcohol-related problems phenotypes (e.g. Alcohol Use Disorders Identification Test (AUDIT) problem score) (rGs: 0.132-0.533, all FDR adjusted p < 0.05). However, the genetic correlations between PTSD phenotypes and alcohol consumption -related phenotypes (e.g. drinks per week) were negatively associated or non-significant (rGs: -0.417 to -0.042, FDR adjusted p: <0.05-NS). For AAs, the direction of correlations was sometimes consistent and sometimes inconsistent with that in EAs, and the ranges were larger (rGs for alcohol-related problems: -0.275 to 0.266, FDR adjusted p: NS, alcohol consumption-related: 0.145-0.699, FDR adjusted p: NS). CONCLUSIONS: These findings illustrate that the genetic associations between consumption and problem alcohol phenotypes and PTSD differ in both strength and direction. Thus, the genetic factors that may lead someone to develop PTSD and high levels of alcohol consumption are not the same as those that lead someone to develop PTSD and alcohol-related problems. Discussion around needing improved methods to better estimate heritabilities and genetic correlations in diverse and admixed ancestry samples is provided.

Examining interactions between polygenic scores and interpersonal trauma exposure on alcohol consumption and use disorder in an ancestrally diverse college cohort.

Introduction: Genetic factors impact alcohol consumption and use disorder (AUD), with large-scale genome-wide association studies (GWAS) identifying numerous associated variants. Aggregate genetic methods in combination with important environmental factors (e.g., interpersonal trauma [IPT]) can be applied to expand our understanding of the ways by which genetic and environmental variables work together to influence alcohol consumption and disordered use. The present study aimed to detail the relationships between genome-wide polygenic scores (PGS) for alcohol phenotypes (i.e., alcohol consumption and AUD status) and IPT exposure as well as the interaction between them across ancestry. Methods: Data were drawn from the Spit for Science (S4S) study, a US college student population, where participants reported on IPT exposure prior to college and alcohol consumption and problems during college (N = 9,006; ancestry: 21.3% African [AFR], 12.5% Admixed Americas [AMR], 9.6% East Asian [EAS], 48.1% European [EUR], 8.6% South Asian [SAS]). Two trans-ancestry PGS were constructed, one for alcohol consumption and another for AUD, using large-scale GWAS summary statistics from multiple ancestries weighted using PRS-CSx. Regression models were applied to test for the presence of associations between alcohol-PGS and IPT main and interaction effects. Results: In the meta-analysis across ancestry groups, IPT exposure and PGS were significantly associated with alcohol consumption (ßIPT = 0.31, P IPT = 0.0002; ßPGS = 0.09, P PGS = 0.004) and AUD (ORIPT = 1.12, P IPT = 3.5 × 10-8; ORPGS = 1.02, P PGS = 0.002). No statistically significant interactions were detected between IPT and sex nor between IPT and PGS. When inspecting ancestry specific results, the alcohol consumption-PGS and AUD-PGS were only statistically significant in the EUR ancestry group (ßPGS = 0.09, P PGS = 0.04; ORPGS = 1.02, P PGS = 0.022, respectively). Discussion: IPT exposure prior to college was strongly associated with alcohol outcomes in this college-age sample, which could be used as a preventative measure to identify students at high risk for problematic alcohol use. Additionally, results add to developing evidence of polygenic score association in meta-analyzed samples, highlighting the importance of continued efforts to increase ancestral representation in genetic studies and inclusive analytic approaches to increase the generalizability of results from genetic association studies.

From genetics to systems biology of stress-related mental disorders.

Many individuals will be exposed to some form of traumatic stress in their lifetime which, in turn, increases the likelihood of developing stress-related disorders such as post-traumatic stress disorder (PTSD), major depressive disorder (MDD) and anxiety disorders (ANX). The development of these disorders is also influenced by genetics and have heritability estimates ranging between ∼30 and 70%. In this review, we provide an overview of the findings of genome-wide association studies for PTSD, depression and ANX, and we observe a clear genetic overlap between these three diagnostic categories. We go on to highlight the results from transcriptomic and epigenomic studies, and, given the multifactorial nature of stress-related disorders, we provide an overview of the gene-environment studies that have been conducted to date. Finally, we discuss systems biology approaches that are now seeing wider utility in determining a more holistic view of these complex disorders.

Mineralocorticoid receptors dampen glucocorticoid receptor sensitivity to stress via regulation of FKBP5.

Responding to different dynamic levels of stress is critical for mammalian survival. Disruption of mineralocorticoid receptor (MR) and glucocorticoid receptor (GR) signaling is proposed to underlie hypothalamic-pituitary-adrenal (HPA) axis dysregulation observed in stress-related psychiatric disorders. In this study, we show that FK506-binding protein 51 (FKBP5) plays a critical role in fine-tuning MR:GR balance in the hippocampus. Biotinylated-oligonucleotide immunoprecipitation in primary hippocampal neurons reveals that MR binding, rather than GR binding, to the Fkbp5 gene regulates FKBP5 expression during baseline activity of glucocorticoids. Notably, FKBP5 and MR exhibit similar hippocampal expression patterns in mice and humans, which are distinct from that of the GR. Pharmacological inhibition and region- and cell type-specific receptor deletion in mice further demonstrate that lack of MR decreases hippocampal Fkbp5 levels and dampens the stress-induced increase in glucocorticoid levels. Overall, our findings demonstrate that MR-dependent changes in baseline Fkbp5 expression modify GR sensitivity to glucocorticoids, providing insight into mechanisms of stress homeostasis.

Increasing the resolution and precision of psychiatric genome-wide association studies by re-imputing summary statistics using a large, diverse reference panel.

Genotype imputation across populations of mixed ancestry is critical for optimal discovery in large-scale genome-wide association studies (GWAS). Methods for direct imputation of GWAS summary-statistics were previously shown to be practically as accurate as summary statistics produced after raw genotype imputation, while incurring orders of magnitude lower computational burden. Given that direct imputation needs a precise estimation of linkage-disequilibrium (LD) and that most of the methods using a small reference panel for example, ~2,500-subject coming from the 1000 Genome-Project, there is a great need for much larger and more diverse reference panels. To accurately estimate the LD needed for an exhaustive analysis of any cosmopolitan cohort, we developed DISTMIX2. DISTMIX2: (a) uses a much larger and more diverse reference panel compared to traditional reference panels, and (b) can estimate weights of ethnic-mixture based solely on Z-scores, when allele frequencies are not available. We applied DISTMIX2 to GWAS summary-statistics from the psychiatric genetic consortium (PGC). DISTMIX2 uncovered signals in numerous new regions, with most of these findings coming from the rarer variants. Rarer variants provide much sharper location for the signals compared with common variants, as the LD for rare variants extends over a lower distance than for common ones. For example, while the original PGC post-traumatic stress disorder GWAS found only 3 marginal signals for common variants, we now uncover a very strong signal for a rare variant in PKN2, a gene associated with neuronal and hippocampal development. Thus, DISTMIX2 provides a robust and fast (re)imputation approach for most psychiatric GWAS-studies.

Intergenerational trauma is associated with expression alterations in glucocorticoid- and immune-related genes.

Offspring of trauma survivors are more likely to develop PTSD, mood, and anxiety disorders and demonstrate endocrine and molecular alterations compared to controls. This study reports the association between parental Holocaust exposure and genome-wide gene expression in peripheral blood mononuclear cells (PBMC) from 77 Holocaust survivor offspring and 15 comparison subjects. Forty-two differentially expressed genes (DEGs) were identified in association with parental Holocaust exposure (FDR-adjusted p < 0.05); most of these genes were downregulated and co-expressed in a gene network related to immune cell functions. When both parental Holocaust exposure and maternal age at Holocaust exposure shared DEGs, fold changes were in the opposite direction. Similarly, fold changes of shared DEGs associated with maternal PTSD and paternal PTSD were in opposite directions, while fold changes of shared DEGs associated with both maternal and paternal Holocaust exposure or associated with both maternal and paternal age at Holocaust exposure were in the same direction. Moreover, the DEGs associated with parental Holocaust exposure were enriched for glucocorticoid-regulated genes and immune pathways with some of these genes mediating the effects of parental Holocaust exposure on C-reactive protein. The top gene across all analyses was MMP8, encoding the matrix metalloproteinase 8, which is a regulator of innate immunity. To conclude, this study identified a set of glucocorticoid and immune-related genes in association with parental Holocaust exposure with differential effects based on parental exposure-related factors.

Genome-wide translational profiling of amygdala Crh-expressing neurons reveals role for CREB in fear extinction learning.

Fear and extinction learning are adaptive processes caused by molecular changes in specific neural circuits. Neurons expressing the corticotropin-releasing hormone gene (Crh) in central amygdala (CeA) are implicated in threat regulation, yet little is known of cell type-specific gene pathways mediating adaptive learning. We translationally profiled the transcriptome of CeA Crh-expressing cells (Crh neurons) after fear conditioning or extinction in mice using translating ribosome affinity purification (TRAP) and RNAseq. Differential gene expression and co-expression network analyses identified diverse networks activated or inhibited by fear vs extinction. Upstream regulator analysis demonstrated that extinction associates with reduced CREB expression, and viral vector-induced increased CREB expression in Crh neurons increased fear expression and inhibited extinction. These findings suggest that CREB, within CeA Crh neurons, may function as a molecular switch that regulates expression of fear and its extinction. Cell-type specific translational analyses may suggest targets useful for understanding and treating stress-related psychiatric illness.

Assessing the Role of Long Noncoding RNA in Nucleus Accumbens in Subjects With Alcohol Dependence.

BACKGROUND: Long noncoding RNA (lncRNA) have been implicated in the etiology of alcohol use. Since lncRNA provide another layer of complexity to the transcriptome, assessing their expression in the brain is the first critical step toward understanding lncRNA functions in alcohol use and addiction. Thus, we sought to profile lncRNA expression in the nucleus accumbens (NAc) in a large postmortem alcohol brain sample. METHODS: LncRNA and protein-coding gene (PCG) expressions in the NAc from 41 subjects with alcohol dependence (AD) and 41 controls were assessed via a regression model. Weighted gene coexpression network analysis was used to identify lncRNA and PCG networks (i.e., modules) significantly correlated with AD. Within the significant modules, key network genes (i.e., hubs) were also identified. The lncRNA and PCG hubs were correlated via Pearson correlations to elucidate the potential biological functions of lncRNA. The lncRNA and PCG hubs were further integrated with GWAS data to identify expression quantitative trait loci (eQTL). RESULTS: At Bonferroni adj. p-value ≤ 0.05, we identified 19 lncRNA and 5 PCG significant modules, which were enriched for neuronal and immune-related processes. In these modules, we further identified 86 and 315 PCG and lncRNA hubs, respectively. At false discovery rate (FDR) of 10%, the correlation analyses between the lncRNA and PCG hubs revealed 3,125 positive and 1,860 negative correlations. Integration of hubs with genotype data identified 243 eQTLs affecting the expression of 39 and 204 PCG and lncRNA hubs, respectively. CONCLUSIONS: Our study identified lncRNA and gene networks significantly associated with AD in the NAc, coordinated lncRNA and mRNA coexpression changes, highlighting potentially regulatory functions for the lncRNA, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.

Author Correction: Genome-wide translational profiling of amygdala Crh-expressing neurons reveals role for CREB in fear extinction learning.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

TWAS pathway method greatly enhances the number of leads for uncovering the molecular underpinnings of psychiatric disorders.

Genetic signal detection in genome-wide association studies (GWAS) is enhanced by pooling small signals from multiple Single Nucleotide Polymorphism (SNP), for example, across genes and pathways. Because genes are believed to influence traits via gene expression, it is of interest to combine information from expression Quantitative Trait Loci (eQTLs) in a gene or genes in the same pathway. Such methods, widely referred to as transcriptomic wide association studies (TWAS), already exist for gene analysis. Due to the possibility of eliminating most of the confounding effects of linkage disequilibrium (LD) from TWAS gene statistics, pathway TWAS methods would be very useful in uncovering the true molecular basis of psychiatric disorders. However, such methods are not yet available for arbitrarily large pathways/gene sets. This is possibly due to the quadratic (as a function of the number of SNPs) computational burden for computing LD across large chromosomal regions. To overcome this obstacle, we propose JEPEGMIX2-P, a novel TWAS pathway method that (a) has a linear computational burden, (b) uses a large and diverse reference panel (33 K subjects), (c) is competitive (adjusts for background enrichment in gene TWAS statistics), and (d) is applicable as-is to ethnically mixed-cohorts. To underline its potential for increasing the power to uncover genetic signals over the commonly used nontranscriptomics methods, for example, MAGMA, we applied JEPEGMIX2-P to summary statistics of most large meta-analyses from Psychiatric Genetics Consortium (PGC). While our work is just the very first step toward clinical translation of psychiatric disorders, PGC anorexia results suggest a possible avenue for treatment.

Analysis of Genetically Regulated Gene Expression Identifies a Prefrontal PTSD Gene, SNRNP35, Specific to Military Cohorts.

To reveal post-traumatic stress disorder (PTSD) genetic risk influences on tissue-specific gene expression, we use brain and non-brain transcriptomic imputation. We impute genetically regulated gene expression (GReX) in 29,539 PTSD cases and 166,145 controls from 70 ancestry-specific cohorts and identify 18 significant GReX-PTSD associations corresponding to specific tissue-gene pairs. The results suggest substantial genetic heterogeneity based on ancestry, cohort type (military versus civilian), and sex. Two study-wide significant PTSD associations are identified in European and military European cohorts; ZNF140 is predicted to be upregulated in whole blood, and SNRNP35 is predicted to be downregulated in dorsolateral prefrontal cortex, respectively. In peripheral leukocytes from 175 marines, the observed PTSD differential gene expression correlates with the predicted differences for these individuals, and deployment stress produces glucocorticoid-regulated expression changes that include downregulation of both ZNF140 and SNRNP35. SNRNP35 knockdown in cells validates its functional role in U12-intron splicing. Finally, exogenous glucocorticoids in mice downregulate prefrontal Snrnp35 expression.

The impact on estimations of genetic correlations by the use of super-normal, unscreened, and family-history screened controls in genome wide case-control studies.

Traditionally, in normal case-control studies of disorder A, the controls are defined as those not developing the disorder. However, in genome wide association (GWA) studies, controls are sometimes (a) unscreened or (b) screened for both disorder A and disorder B, producing super-normal controls. Using simulations, we examine how the observed genetic correlations between two disorders (A and B) are influenced by the use of unscreened, normal, and super-normal controls. Normal controls produce unbiased estimates of the genetic correlation. However, unscreened and super-normal controls both bias upward the genetic correlations. The strength of the bias increases with increasing population prevalences for the two disorders. With super-normal controls, the absolute magnitude of bias is stronger when the true genetic correlation is low. The opposite is seen with the use of unscreened controls. Adding screening of first-degree relatives of controls substantially increases the bias in genetic correlations with super-normal controls but has minimal impact when controls are screened only for the relevant disease.

Genome-wide association study of shared liability to anxiety disorders in Army STARRS.

Anxiety disorders (ANX), namely generalized anxiety, panic disorder, and phobias, are common, etiologically complex syndromes that show increasing prevalence and comorbidity throughout adolescence and beyond. Few genome-wide association studies (GWAS) examining ANX risk have been published and almost exclusively in individuals of European ancestry. In this study, we phenotyped participants from the Army Study To Assess Risk and Resilience in Servicemembers (STARRS) to approximate DSM-based ANX diagnoses. We factor-analyzed those to create a single dimensional anxiety score for each subject. GWAS were conducted using that score within each of three ancestral groups (EUR, AFR, LAT) and then meta-analyzed across ancestries (NTotal = 16,510). We sought to (a) replicate prior ANX GWAS findings in ANGST; (b) determine whether results extended to other ancestry groups; and (c) meta-analyze with ANGST for increased power to identify novel susceptibility loci. No reliable genome-wide significant SNP associations were detected in STARRS. However, SNPs within the CAMKMT gene located in region 2p21 associated with shared ANX risk in ANGST were replicated in EUR soldiers but not other ancestry groups. Combining EUR STARRS and ANGST (N = 28,950) yielded a more robust 2p21 association signal (p = 9.08x10-11 ). Gene-based analyses supported three genes within 2p21 and LBX1 on chromosome 10. More powerful ANX genetic studies will be required to identify further loci.